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ABSTRACT: In their 1997 paper, CDF measured a e ff, the normalization factor that relates 
the cross section for double parton scattering to the product of the inclusive cross sections 
of the two individual scatters, in a model in which they are assumed to be independent. 
In his 2007 paper, Treleani pointed out that CDF used a non-standard definition, in which 
the double parton scattering cross section corresponds to exactly two scatters, rather than 
the more conventional one in which it is the inclusive two-scatter cross section. He also 
estimated the correction from one definition to the other, to give a corrected value of a e g. 
Treleani's form would be correct under the assumption that CDF were able to uniquely 
identify and count the number of scatters in an event, which is certainly not the case. In 
this publication we consider CDF's event definition in more detail to provide an improved 
correction. 
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1 Introduction 

Due to the composite nature of hadrons, it is possible to have multiple parton scatterings, 
i.e. events in which two or more distinct parton interactions occur simultaneously in a 
single hadron-hadron collision. At fixed final state invariant masses, such cross sections 
tend to increase with collision energy because partons with successively lower momentum 
fraction x, hence rapidly increasing fluxes, are probed. Therefore, the question of multiple 
parton interactions (MPI) in a single hadronic collision has rapidly moved from a theoretical 
curiosity, when double parton interactions was for the first time observed by AFS experiment 
at t/s = 63 GeV [1] to a critical issue at the LHC. Multiple parton interactions at the 
LHC give rise to different effects, among others a substantial increase of the unavoidable 
background to most observables used for the search of new physics [2-7]. For this reason, 
multiple parton scattering has taken on considerable importance in recent years, since a 
variety of new and improved Monte Carlo models of underlying event physics rely on it 
[8-17]. Unfortunately the process cannot be estimated in a straightforward way, due to 
the lack of knowledge of the non-pert urbative physics, therefore these models rely on the 
experimental input. Especially, the experimental value of a e ff, the normalization factor 
that relates the cross section for double parton scattering to the product of the inclusive 
cross sections of the two individual scatters, has the potential to act as a strong constraint 
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on models of multiple parton scattering and, in particular, their models of the transverse- 
space distribution of partons in hadrons. The value of a e g is also crucial for calculations 
of double pair production based on the framework of the /cj-factorization approach, see for 
example [18]. 

CDF's measurement of the double parton scattering cross section [19], still one of the 
best available, therefore has considerable significance not only for MPI models but also as a 
standard candle for new [20-22] and planned measurements [23]. CDF's measurement was 
better than any that came before it, because it avoided almost all reliance on a Monte Carlo 
description of their final states and on theoretically-calculated cross sections. Instead it 
made an ingenious direct extraction of a e g by defining an event selection sensitive to double 
parton scattering and comparing the rate of these events from beam crossings with a single 
vertex (assumed to be double parton scattering within one proton-antiproton collision) 
and with two vertices (assumed to be single-parton scatterings within two independent 
proton-antiproton collisions). The main assumption that their extraction relies on is that 
the final state of two scatters in the same proton-antiproton collision is identical to that 
of two scatters in different proton-antiproton collisions. CDF defined the cross section 
normalization factor a e jj through the equation 

&a<?b ... 

Vab;2 = , (1-1) 

eff, CDF 

where a a b- t 2 is the cross section for a colliding proton-antiproton pair to have exactly two 
scatters of types a and b. We choose the notation for exclusive cross sections, where the 
process is described by the subscript containing the type of individual parton processes and 
the number of scatters after a semi-colon. Here, we assume that the two processes a and b 
are different, i.e. distinguishable 1 . a a f> are their inclusive cross sections. To distinguish this 
from the definition more commonly used in theoretical studies, we write it as (J e ff t cDF from 
here on. CDF's final value was 

eff, CDF = (14.5 ± Mt\l) mb. (1.2) 

Since they used the definition of (1.1), which refers to exactly two scatters, they made 
a correction to account for the fact that a fraction of their events (which they estimated 
to be 17^8%) came from triple-parton scattering events. In their analysis CDF used this 
estimation to re-scale the number of the accepted events by the factor of 0.83^q'q4. The 
ratio of these two numbers is equal to the ratio of the exclusive triple and double parton 
scattering cross sections producing the same final state: 



&ab:2 83 



(1.3) 



From the theoretical point of view, it is more convenient to define a e ff through an 
analogous formula, but for the inclusive double-scattering cross section, 

<? a Vb .s 

(Tab = , (1-4) 

O 'eff 



1 CDF use an experimental method to verify that this is the case with their event selection, even though 
there could be some overlap in principle. 



- 2 - 



since, as shown in Section 3, with this definition in the assumption of independence of 
individual scatters, a e jj depends only on properties of the colliding hadrons and not on the 
scattering types or cross sections. 

The non-standard definition used by CDF has been pointed out for the first time by 
Treleani in his publication [24]. Based on the theoretically-pure (parton-level) situation, 
assuming that one could measure, for a given event, whether it came from two scatters of 
definite types a and b or three, one of type a and two of type b, Treleani calculated the 
correction from (J e ff,CDF to a e jj. Using the formula in his paper, one obtains the value 2 



The purpose of this publication is to point out that this extraction is over-simplified - 
the indirectness of CDF's measurement means that they are far from being in this pure 
(parton-level) situation. 

The paper is organised as follows. We begin by recapping the salient points of CDF's 
measurement and event selection in Sect. 2. In Sect. 3 we consider how to calculate the two- 
and three-scatter cross sections with this event selection. We show, in Sect. 4, that further 
input is needed to relate these to a e g. In Sect. 5 we make Monte Carlo estimates of this 
input and its uncertainty. One important point to note already is that we only need ratios 
of closely-related cross sections, so we hope not to be too sensitive to details of the Monte 
Carlo event generation so that the uncertainty is fairly small. Finally, in Sect. 6 we put 
this input together with CDF's measurement to give a final value for a e jj. In appendices 
we give more details backing up our Monte Carlo evaluation of the correction factor and 
its uncertainty. 

2 CDF's experimental measurement 

The intricacies of the clever extraction of a e ff using the single- and double-vertex data 
will not need to concern us here. What will be important is the definition of the final 
state containing a direct photon and exactly three jets. The photon was measured within 
pseudo-rapidity acceptance \rj\ < 0.9 and was required to have Et^ > 16 GeV. Jets were 
reconstructed from all objects in the calorimeter region, \rj\ < 4.2, using a cone algorithm 
CDFJetClu with cone radius 0.7 [25]. Jets were ordered according to their Et so that 
Eti > Et2 > Et3 and were required to be separated in rj — <j) plane from each other by 0.7 
and from the photon by 0.8. Events were accepted if all three jets had Et > 5 GeV and 
simultaneously the second and the third hardest jets had Et < 7 GeV. 

In what we will call the pure parton picture, one would therefore have one of the jets 
also above 16 GeV, equal and opposite to the photon in transverse momentum both coming 
from one scatter, and the other two jets equal and opposite to each other, from the other 
scatter. However, with 5 GeV jets, one is far from this theoretically pure situation. CDF 
spent some time investigating the properties of such low Et jets and concluded that they 
are sufficiently correlated with the underlying parton dynamics to enable their measurement 

2 Which in [24] is written as "a e ff ~ limb". 



10.3 mb. 
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but that there was a great deal of smearing and creation of jets 'from nothing'. We certainly 
cannot therefore rely on this simple parton picture. In fact, CDF estimated that 75% of 
their event sample came from events in which the photon and two of the jets came from 
one scattering and one jet from the other. They stressed throughout that the numerator 
a a a b hi Eq. (1.1) is a shorthand for the sum over all separations into two scatters of the 
source of their three jets. In this spirit, one could write more precisely 

07+ljetC2jets + 0"7+2jetsO"ljet fn 1 \ 

CT 7 +3jets;2 = • (2.1 J 

&eff,CDF 

Note that the contribution containing cr 7 _|_oj e ts, which should in principle also appear here, 
was found to be negligible [19]. 

The data were corrected for the trigger efficiency, but not for any other detector effects, 
but we argue below (see Sec. 5) that, since we will only need ratios of closely-related cross 
sections, they should largely cancel. 

3 Parton level correction 

In order to understand the correction at the parton level from (7 e ff t cDF to a e ff, we write 
down the expression for the cross section cr ab ;2 according to Eq. (2.1) in a simple eikonal 
model. We assume that the dependence of the parton distributions on the two-dimensional 
impact parameter, b, and longitudinal momentum fraction, x, factorize. We make extensive 
use of the overlap function A(b), normalized such that 



d 2 6 ^(6) = 1. (3.1) 

To set the scene, we mention that the cross section for exactly n scatters of a type b is 
given by 



<Jb;n = J d 2 b^a b A(b)) n exp^-a b A(b)y (3.2) 



That is, the cross section is obtained by integrating over all values of impact parameter the 
probability of the scatters: each has a probability a b A(b); there are n of them and they are 
independent, giving the power of n; and they are all of the same type, giving the n! factor. 
Finally, the exponential gives the probability that there are no other scatters (of type b). 
It is straightforward to check that the inclusive cross section is given by 



0"6 



^2<na b;n . (3.3) 



For a total of n scatters: exactly one of type a and n— 1 of type b, with a and b assumed 
distinguishable, we obtain in the same way 

<Jab-,n = f d 2 b[a a A{b)) l _ ; (* b A (b)) expj - (a a + a h ) A(b) \ . (3.4) 

In all that follows, we assume that one of the two scattering types, let us say a, has a small 
cross section, so we can drop a a from the exponent. We therefore have, for the exactly 



- 4 - 



two-scatter cross section 



<?ab;2 = <r a a b J d 2 b(^A(b)j expj -a b A(b) j , 



(3.5) 



and hence 

O eff, CDF = —. — 2 7 S • ( 3 - 6 ) 

Jd 2 b[A(b)) expl-a b A(b)\ 
If, instead, we define a e g through the inclusive two-scatter cross section 

<7ab = ^2( n ~ l )Oab;n = (?a<7b j d 2 b^A(b^j , (3.7) 

we obtain 

J -^2 > (3-8) 
d 2 b(A(b)j 

independent of the cross sections for the individual processes selected. 

It is clear that if the cross section for process b is small, the two definitions, Eqs. (3.6) 
and (3.8), become identical. Keeping the first correction to this, we obtain 

o eff, cdf ~ o e g + TZa b , (3.9) 

where 



d 2 b 



n=^ ^ (3.io) 



/ ->\ 2 



.PMi(!))J 

is a function only of the shape, but not the size, of the overlap function. For example, for 
a Gaussian matter distribution we have 7Z = 1.33, for a matter distribution based on the 
electromagnetic form factor, as used in Refs.[10, 15], we have 1Z = 1.46, for an exponential 
matter distribution, we have 7Z = 1.78 and for a 'black disk' we have 7Z = 1.26 (we mention 
these numbers for illustration, but they are not needed for our extraction). 
Treleani's idea is to obtain TZa b term from the rate of 3-scatter events, 

<Tab;3 = J d 2 b(a a A{b)) ^ (a b A(b)) 2 exp|-a 6 A(b)\ . (3.11) 

With the same accuracy as we have just used to obtain Eq. (3.9), we can replace the 
exponential by unity and obtain 

a ab ;3 « \faol j d 2 b[A(b)f « ^a ab , 2 ^ . (3.12) 

By insertion of 7Z from (3.12) into (3.9), we completely reproduce the correction suggested 
in [24]: 

°eff,CDF ~ <? e ff f 1 + 2°^- ) . (3.13) 
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4 Jet level correction 



We turn now to a calculation of the same quantities, with the same accuracy, but at the jet 
level using the event definition of CDF. By jet level, we mean that we take account of the 
subsequent evolution of the parton-level event into a multi-parton system and thence into 
the hadrons that are reconstructed as jets in the detector. These steps smear the results 
considerably, due to the emission of extra jets into the event, the recoil of the primary 
jets from them, the smearing of jet momenta by hadronization effects, the loss of hadronic 
energy out of the jet, the merging of nearby jets, etc. For jet physics with a threshold of 
5 GeV, all of these effects are large. Not only are the jets smeared, but it becomes impossible 
to match them up and know which jets came from which scatters. 

Returning to Eq. (2.1), i.e. including the fact that the observed three jets can come 
from two scatters in one of two ways, we obtain 

0- 7 +3jets;2 = / d ^ (f-y+ljet A(b)j (a 2 jeta A(b)j 

+ (tf 7 +2jcts A(bj^ (ffijct A(bj) exp|-a jct s^(&)| • (4.1) 

Note that the cross section in the exponent is a third type of scattering: in both separations 
of events, we veto all scatters with any jets above 5 GeV. Therefore Ojets is the cross section 
for a scatter to produce at least one jet with Ej- > 5 GeV. With the same accuracy as the 
argumentation made above, we therefore have 

&eff,CDF ~ O- e jj + TZa-^ts- (4.2) 
On the other hand, the three-scatter contribution to this cross section is given by 

o- 7 +3jets;3 = ^ J d 2 b (cr 7+ ij e t A(b)^j (aij et A(bfj (o-ijct .4(6)) expj-o- jets A(b)\ . (4.3) 

That is, we continue to neglect the 7 + Ojets term, so three jets from three scatters can only 
come about from each of the extra scatters contributing exactly one jet. Approximating 
this in the same way as above and putting everything together we obtain 



, . o q- 7 +3jets;3 ^ 7 +ljet<72jets + <7 7+ 2jets^ljet <Tjets \ 

&eff,CDF ~ 0- e ff[ 1 + Z X X (4.4 J 

<7 7 +3jets;2 0" 7 +ljet<7ljet °"ljet/ 

a I l ! Q ^T+SjctsjS / gjjets | *7 7 +2jcts \ (Tjcts \ ^ ^ 

0"7+3jcts;2 \ CTljet 07+ljet / Cljct/ 

! 1 • 2 a ~- :U ^.l). (4.6) 



0"7+3jets;2 

where in the last equation we define the correction factor 



j. _ / <72jets + 0'7+2jet S \ gjets ^ ^ 



°~ljet 07+ljet J 0"ij e t 

Note that at the parton level, in which scatters can be identified and counted perfectly and 
jets are replaced by partons, the correction factor / is equal to unity since each of the two 
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new ratios that have appeared in Eq. (4.4) is equal to one, since in this case <7ij e t, 0"2jets 
and Ojets are all equal and <7 7 +2j e ts i s negligible. In general, however, the first new factor 
can be expected to be smaller than unity and the second one larger than unity so it is not 
clear the direction of the overall effect. 

We do not believe it is possible to extract these ratios from numbers in the CDF paper 
alone. Nevertheless, since they are ratios of closely-related cross sections, we may hope 
that they are considerably better predicted than the cross sections themselves. In the next 
section we use Monte Carlo data to extract the correction factor and make an estimate of 
its uncertainty. 

5 Estimate of the jet correction factor / 

Before proceeding to the numerical estimate, we comment on the accuracy required. Since 
the ratio (1.3) has an uncertainty of ~ ^50%, an uncertainty of the correction factor / that 
multiplies it of order 20% or less would be ample, leaving the result for the effective cross 
section dominated by the experimental uncertainties. It would be unreasonable to aim for 
a significantly higher accuracy than this, since 05(5 GeV) ~ 0.2 and this analysis relies on 
the leading order Monte Carlo generators. 

In order to calculate the correction factor / and more reliably estimate its accuracy, we 
used three different Monte Carlo event generators Herwig++ 2.5.2 [26, 27], fHerwig 6.510 [28, 
29] and Pythia 6.4.26 [30] to produce events of both 7 + jets and pure QCD jets types. Each 
generated event was handed over to the Rivet package [31] to be analyzed. This ensured 
that the computation of observables is exactly the same for each generator. Since the CDF 
analysis was not available among the standard Rivet's analyses we have implemented its 
event selection criteria (see Section 2) into the package. For jet finding, we used the default 
settings of CDFJetClu jet algorithm as implemented in FastJet 2.4.2 [32, 33] with a cone 
radius R = 0.7, which according to [25] was used in the original CDF analysis. 

In the CDF analysis, jet cuts depend on the order of jet in the event with respect to 
its transverse energy. The fact that the leading jet could be above or below the 7 GeV, 
while the two trailing jets have to be below this threshold, induces an additional correlation 
between the two scatters, in principle. The sum over the two divisions of the jet origin 
between the two scatters should be extended to include a sum over all assignments of the 
leading jet between the two scatters, and whether it is above or below 7 GeV. However, the 
part of the cross section coming from events in which the leading jet comes from the QCD 
scattering and not from the photon production was found to be tiny, about 1.5%. Moreover, 
it is about the same fraction in the numerator and the denominator for the appropriate 
ratio, so neglecting these events really has a negligible effect. Therefore, the highest Ey jet 
in 7 + jets is required to be above 5 GeV and all other jets to be between 5 GeV and 7 GeV. 

Since we want to make predictions for the final state of a single scattering, we switch 
off all MPI effects, but leave all other generator options on and at their default values. As 
we have emphasised, jets of such low transverse momentum are strongly smeared and only 
weakly correlated with the partonic scattering that produced them. Correspondingly, there 
is a significant probability that jets that enter our acceptance may be initiated by partonic 
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scatterings with transverse momentum well below that acceptance. In order to estimate the 
amount of this migration, which also gives an idea of the contamination of the jet sample by 
very soft physics, we divide the calculated cross sections into two separate parts, which we 
call Soft and Hard, according to the matrix element pt- For the jets (7 + jets, respectively) 
production cross sections the Hard part is defined with pt > 2 GeV (pt > 10 GeV) and the 
Soft 3 with 0.5 < p t < 2 GeV (5.0 < p t < 10 GeV). 

The obtained cross sections, their ratios and values of the correction factor for each 
generator are presented in Table l 4 . The first observation is that the 7 + jets sample for 
Et~{ > 16 GeV is well behaved: there is almost no migration from hard processes below 
10 GeV. However, in the case of the QCD scattering cross sections we see a significant 
amount of migration from below 2 GeV. We conclude that these events are somewhat less 
likely to be well modeled, which is also reflected by significantly different cross sections 
obtained from different generators. Nevertheless, the amount of migration is similar in each 
of these cross sections, therefore the mis-modeling cancels to some extent in their ratios 
and hence the correction factor is reasonably predicted. This is particularly evident when 
comparing the correction factors obtained from Herwig++ and Pythia 6. Despite significant 
differences in the cross sections, the values of the correction factors are very close. The 
reason why this coefficient is different in the case of fHerwig will be explained later in this 
section. 
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a [mb] 


Hard 


Soft 


Sum 
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Sum 


Hard 
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9.16 




3.16 
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5.33 


6.61 


11.94 
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2.51 


9.44 


&2jets 
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0.77 
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0.70 
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0.07 
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0.103 






0.076 
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1.426 






1.426 






1.383 




&~f + 2je.ts 






0.258 






0.300 






0.246 




/ 


0.458 


0.575 


0.445 


favg. 


0.493 



Table 1. The calculated cross sections, their ratios and values of the correction factor for the 
default settings of each generator Herwig++, fHerwig, and Pythia 6. 



We have already begun to address the question of how the results depend on the 
Monte Carlo modeling using three different event generators with their default settings. 
In addition, we use the Herwig++ generator with the default settings altered in order to 
determine how the details of the simulation affect the results. More details are given in the 

3 It is worth noting that default setting of Pythia 6 does not produce any events with pt < 1.0 GeV. 
4 The statistical errors are negligible, therefore we suppress them in the Table. 
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appendices. The most important effects for our studies are the order of as (1-loop and 2- 
loops) used in the simulation, width of the Gaussian distribution of the intrinsic transverse 
momentum kr of the interacting partons (we studied three values 0, 1 and 2 GeV) and 
the parton distribution functions. In addition, since CDF did its studies at calorimeter 
level and we do it at particle level, we believe that there is value in using a different jet 
algorithm (CDFJetClu, PxCone[34] 5 and Anti-k^[35] as implemented in FastJet 2.4.2), to 
see how dependent we are on these fine details. In Table 2 we present results obtained 
using three different PDF sets (MRST98 [36], CTEQ6L1 [37] and MRST LO** [38]) and 
two different orders of as- We can see that the impact of the PDF on the result is small. 
Similarly, the order of as has little effect on the / factor. The jet clustering algorithms, see 



PDF 


MRST98 


CTEQ6L1 


MRST LO** 


a s 


1-loop 


2-loops 


/ 


0.477 


0.447 


0.458 


f 


0.476 


0.458 



Table 2. The correction factors obtained using Herwig+- 1- with three different PDF sets MRST98, 
CTEQ6L1 and MRST LO** (default in Herwig+- H) and two different orders of as, 1-loop and 2-loops 
(default in Herwig+- H). 

Table 3, have slightly bigger influence on the outcome but still smaller then the uncertainty 
coming from the different Monte Carlo models, see Table 1. 



Jet algorithm 


CDFJetClu 


PxCone 


Anti-kr 


/ 


0.458 


0.512 


0.525 



Table 3. The correction factors obtained using Herwig+- 1- with three different jet algorithms 
PxCone, Anti-ky and CDFJetClu (used in the CDF analysis). 

By far, the dominant effect is due to the intrinsic kr modeling, therefore we have 
studied its influence in more detail using all three generators. The results from Table 4 
explain why the results obtained using default settings of Pythia 6 and Herwig++ are very 
similar (see Table 1). This is because the intrinsic momentum in both generators was tuned 
to experimental data and have by default similar value kr ~ 2 GeV, while in fHerwig it 
was not tuned to the data and by default is equal to GeV. Therefore, in this respect 
results from Herwig++ and Pythia 6 should be trusted more then from fHerwig. We also 
see (Table 4) that all generators provide a similar value of / for the same kr value. 



/ 


far = 0.0 GeV 


far = 1.0 GeV 


k T = 2.0 GeV 


Herwig++ 


0.648 


0.582 


0.465 


fHerwig 


0.575 


0.619 


0.564 


Pythia 6 


0.620 


0.590 


0.445 



Table 4. The correction factors obtained using three generators and three different values of 
intrinsic far. 

5 The minimal transverse energy threshold was changed from the default value of 0.5 GeV to 0.1 GeV. 
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The more detailed results including the cross sections and their ratios are included in 
Appendix A. 

As an estimate of the correction factor we take the average value of / obtained from the 
three different event generators with their default settings (see Table 1). The systematic 
error of the estimation is taken as a half of the difference between the maximum and 
minimum value of / caused by the effects studied in this section. Therefore, the final result 
is 

favg. = 0.49 ±0.10 . (5.1) 

Which indeed, as we anticipated from the estimation of a s (5 GeV), is around 20% of the 
correction factor. 

As a final cross-check, we calculate the fraction of 7 + 3jets events that come from 
a 7 + 2jets collision plus a ljet collision for which the experimental value was quoted as 
« 75%[19]. We obtain 80% in Herwig++, 74% in fHerwig and 76% in Pythia 6. Taking into 
account the inherent uncertainties in jet physics at 5 GeV, the fact that we are working at 
particle level and they work at uncorrected detector level, and the accuracy we are aiming 
for in the final correction factor, we consider this to be very good agreement with the 
experimental number. 



6 Result 

Using the result of the previous section and the numbers in CDF's paper, we are ready 
to extract a value for a e ff. However, we first make one final comment, concerning the 
systematic errors. In CDF's analysis, the correction for triple-scattering events makes one 
of the biggest single contributions to the final systematic error, since they subtract the 
roughly estimated number of triple-scattering events off the number of double-scattering 
events. We revise CDF's method and define a new effective cross section, (T e ff,CDF, which 
is the value of a e g they would have obtained by keeping triple-scattering events in their 
sample and thus we avoid this source of the systematic error. Their master formula is 



Ndi (A dp \ 

U ^ aDF= N^ P {A^) Rc(7NSD > (fU) 



with 



N D i = 1060 ± 110 ± 110, (6.2) 

N DP+TP = 8865 ± 430 ± 150, (6.3) 

R DP = 0.83 ± 0±g;g|, (6.4) 

Ndp = N DP+TP x R DP = 7360 ± 3601™ (6.5) 

A D p/Adi = 0.958 ± ± 0, (6.6) 

R c = 2.06 ±0.02+°;^, (6.7) 

a NSD = (50.9 ±0± 1.5) mb. (6.8) 

where we have consistently written the first error as statistical and the second as system- 
atic, even when they are assumed to be zero. Ndi stands for the number of double hadron 
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interactions (DI) identified using the vertex detector, Np>p is the number of pure exclu- 
sive double parton scattering (DP) events, whose fraction within all measured multiple- 
scattering events {Ndp+tp) is Rdp- Factors Ap>p and Ap>i characterize the acceptances 
of the appropriate kinematic selections, except the requirements connected to vertex re- 
construction. Coefficient R c represents the ratio between the number of beam crossings in 
which the detector found only one vertex and the number of beam crossings in which the 
detector reconstructed exactly two vertices. All beam crossings are taken into consideration 
in which any kind of non-single-diffr active (NSD) inelastic proton-antiproton collision was 
detected. The appropriate cross section (Jnsd is also given above. 

Our new effective cross section has the same definition as in (6.1) but with Np>p replaced 
by N DP+TP , 

®eff,cdf = A Nm ( 4^ J RcVNSD = (12.0 ± 1.4+}|) mb. (6.9) 
1\£)P+TP V A DI J 

Note that the fractional systematic errors are indeed significantly smaller. 
Finally, we are ready to calculate a e g. In terms of <J e ff,CDF it is given by 



Given that correction factor f avg , is 0.49 ± 0.10, the factor in square brackets turns out to be 
very close to zero indicating that the difference between the inclusive measurement and the 
true inclusive cross section for double parton scattering is very small. The uncertainty on 
the triple-scattering event fraction O.lT^Qgg can be neglected with respect to the uncertainty 
of the 2f avg . — 1 term. Our final result for the effective cross section is 



We see that despite the additional uncertainty coming from the correction factor, the sys- 
tematic uncertainty is smaller and more symmetrical than on the CDF's final value. 

7 Conclusions 

The CDF measurement of double parton scattering [19] used a non-standard definition of 
a e ff that makes this important quantity process dependent. Therefore, the value provided 
by the experiment <J e ff,CDF = (14.5 ± 1.7^2'g) m b is not suitable for comparisons with other 
measurements or as input for theoretical calculations or Monte Carlo models. The non- 
standard definition used by CDF has been pointed out and corrected for the first time by 
Treleani in his publication [24]. Based on the theoretically-pure (parton-level) situation, [24] 
estimated the inclusive (process independent) a e jf = 10.3 mb. This result would be correct 
under the assumption that CDF were able to uniquely identify and count the number of 
scatters in an event, which is certainly not the case. In this publication we have considered 
CDF's event definition in more detail to provide an improved correction leading to 




(6.10) 



a eff = (12.0 ± lAt\i) mb. 



(6.11) 



a eff = (12.0±1.4^J) mb. 



(7.1) 
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It is worth noting that both statistical and systematic uncertainties have decreased, since the 
additional uncertainty of our correction factor is much smaller than the avoided uncertainty 
stemming from the triple scattering removal done originally by CDF. 

The obtained value of a e g serves as a constraint on the Monte Carlo models since the 
recent tunes of MPI models to the LHC data predict its value to be between 25 — 42 mb [39]. 
This inconsistency between theory and experiment can also be seen from Eq. 4.2, if one 
uses the a^ts around 15 mb, see Table 1, and 1Z from Section 3. This behavior suggests 
that the overlap function used in the MC models is oversimplified and should be improved, 
for example, by including ^-dependence [40-42]. This value can also help to understand 
the recent results from the LHCb experiment [20] (see page 23, Fig 10). The experimental 
results for a e g extracted from the production of J /if; mesons together with an associated 
open charm hadron and from double open charm hadron production are different by a factor 
of between two and three. 
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A Detailed Results 

A.l Intrinsic kr dependence 





Herwig++ 


fHerwig 


Pyth 


ia 


6 


a [mb] 


Hard 


Soft 


Hard 


Soft 
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Soft 


&ljet 


5.13 




1.40 


5.33 




6.61 


4.79 




0.06 


^2 jets 


0.65 




0.26 


0.54 




0.70 


0.66 




0.00 


CT jets 


8.67 




2.05 


8.72 




8.31 


8.04 




0.06 


a [nb] 








Cy+ljet 


5.38 




0.06 


3.41 




0.16 


4.46 




0.08 


cr 7+2 jets 


1.39 




0.01 


1.02 




0.04 


0.92 




0.15 


&2jets 
&ljet 


0.139 




0.103 




0.136 




O jets 
&ljet 


1.641 




1.426 




1.668 




CT 7+2 jets 
V-t+ljet 


0.256 




0.300 




0.236 




f 


0.648 


0.575 


0.620 



Table 5. The calculated cross sections, their ratios and the final correction factors for the intrinsic 
Ut RMS = 0.0 GeV for three MC generators. CDFJetClu jet algorithm was used. 
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fHerwig 
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ia 


6 


a [mb] 


Hard 


Soft 


Hard 


Soft 


Hard 


Soft 


&ljet 


5.72 




0.86 


5.60 


6.96 


4.78 




0.10 


C 2jets 


0.64 




0.08 


0.54 


0.72 


0.75 




0.00 


&jets 


9.37 




1.11 


9.01 


8.88 


8.01 




0.10 


a [nb] 








Cj+ljet 


5.50 




0.09 


3.46 


0.11 


4.48 




0.08 


cr 7+2 jets 


1.43 




0.01 


1.13 


0.06 


0.92 




0.07 


""2 jets 

°\jet 
° jets 

V\jet 
C r j+2jets 


0.109 
1.591 
0.257 


0.100 
1.424 
0.335 


0.136 
1.668 
0.218 


f 


0.582 


0.619 


0.590 



Table 6. The calculated cross sections, their ratios and the final correction factors for the intrinsic 
Ut RMS = 1.0 GeV for three MC generators. CDFJetClu jet algorithm was used. 
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Herwig++ 


fHerwig 


Pyth 


ia 6 


f — ul — 

(7 [mb] 


Hard 


Soft 
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Hard 


Soft 


a ljet 


9.73 


3.79 


6.23 


10.74 


6.93 


2.51 


C 2jets 


0.65 


0.13 


0.72 


0.81 


0.72 


0.00 


Gjets 


14.73 


4.38 


10.87 


12.99 


10.54 


2.52 


a [nb] 








Cy+ljet 


5.51 


0.05 


3.43 


0.17 


4.47 


0.08 


cr 7+2 jets 


1.48 


0.03 


1.08 


0.05 


1.05 


0.07 


""2 jets 
<J\jet 
O jets 

G\jet 
&"l + 2jets 

G-l+ljet 


0.058 


0.090 


0.076 


1.414 
0.271 


1.404 
0.312 


1.383 
0.246 


f 


0.465 


0.564 


0.445 



Table 7. The calculated cross sections, their ratios and the final correction factors for the intrinsic 
Ut RMS = 2.0 GeV for three MC generators. CDFJetClu jet algorithm was used. 



A. 2 Jet algorithm dependence 
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3.91 


8.89 




3.21 


C 2jets 


0.62 




0.15 


0.88 




0.24 


0.70 
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@jets 


13.87 




3.70 


16.02 




4.56 


13.38 




3.98 


a [nb] 








a J+ljet 


5.66 




0.03 


5.28 




0.04 


2.89 




0.01 


cr 7+2 jets 


1.46 




0.01 


1.53 




0.02 


0.83 




0.00 


^2 jets 
<?ljet 


0.063 




0.076 




0.078 




O jets 
Cljet 


1.426 




1.390 




1.434 




&~f-\-1jets 
V-t+ljet 


0.258 




0.292 




0.288 




f 


0.458 


0.512 


0.525 



Table 8. The calculated cross sections, their ratios and the final correction factors in dependence 
on the jet clustering algorithm using Herwig++ generator. 
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CDFJetClu 


PxC 


one 
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f — ul — 

(7 [mb] 


Hard 


Soft 


Hard 


Soft 


Hard 


Soft 


a ljet 


5.33 


6.61 


6.18 


8.30 


5.73 


9.58 


C 2jets 


0.54 


0.70 


0.78 


1.25 


0.71 


1.96 


&jets 


8.72 


8.31 


9.85 


10.72 


9.27 


13.84 


a [nb] 








Cy+ljet 


3.41 


0.16 


3.07 


0.11 


1.64 


0.02 


cr 7+2 jets 


1.02 


0.05 


1.05 


0.06 


0.45 


0.03 


""2 jets 
<J\jet 
O jets 

G\jet 
&"l + 2jets 

G-l+ljet 


0.103 


0.140 


0.175 


1.426 
0.300 


1.421 
0.348 


1.510 
0.291 


f 


0.575 


0.693 


0.704 



Table 9. The calculated cross sections, their ratios and the final correction factors in dependence 
on the jet clustering algorithm using fHerwig generator. 
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&ljet 


9.62 
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11.41 




3.98 


8.84 




2.50 


(^2 jets 


1.17 




0.00 


1.51 




0.00 


1.16 




0.00 


@jets 


15.04 




2.92 


17.20 




4.02 


13.81 




2.55 


a [nb] 








&j+ljet 


4.47 
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4.21 




0.07 


1.90 




0.02 


C7+2 jets 


1.05 




0.07 
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0.07 


0.45 




0.07 


0~2jets 
V\jet 
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0.085 




0.075 




® jets 
0\jet 


1.383 




1.338 




1.378 




C r j + 2jets 
V~f + ljet 


0.246 




0.258 




0.270 




f 


0.445 


0.459 


0.475 



Table 10. The calculated cross sections, their ratios and the final correction factors in dependence 
on the jet clustering algorithm using Pythia 6 generator. 
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A. 3 PDF dependence 
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&jets 
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3.70 
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C2jets 
&ljet 
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O jets 
&ljet 
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°~, + ljet 
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0.258 


f 


0.4 


77 


0.447 


0.4 


58 



Table 11. The calculated cross sections, their ratios and the final correction factors in dependence 
on the parton distribution function used in Herwig++ generator. 
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0.00 


1.17 
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&jets 


15.04 


2.92 


15.04 


2.92 


a [nb] 








4.47 


0.08 


6.32 


0.09 




1.05 


0.07 


1.29 


0.19 


Oljets 
<y\jet 
O jets 
<J\jet 
<?~( + 2jets 
<?~/ + ljet 


0.076 
1.383 
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1.436 
0.232 


f 


0.445 


0.468 



Table 12. The calculated cross sections, their ratios and the final correction factors in dependence 
on the parton distribution function used in Pythia 6 generator 
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A.4 Order of as dependence 
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5.23 
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&2jets 
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0.15 


&jets 


20.96 
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Vy+ljet 


6.53 


0.05 


5.66 


0.03 


(J 7+2jets 


1.77 


0.04 


1.46 


0.01 


f ljet 


0.062 
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Gljet 


1.412 


1.426 


&~1 + 2jets 
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0.275 


0.258 
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0.476 


0.458 



Table 13. The calculated cross sections, their ratios and the final correction factors in dependence 
on the order of as used in Herwig++ generator. 
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